Magneto-Josephson effects and Majorana bound states in quantum wires 
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A prominent signature of Majorana bound states is the exotic Josephson effects they produce, 
the classic example being a fractional Josephson current with An periodicity in the phase differ- 
ence across the junction. Recent work established that topological insulator edges support a novel 
'magneto- Josephson effect', whereby a dissipationless current exhibits 47r-periodic dependence also 
on the relative orientation of the Zeeman fields in the two banks of the junction. Here, we explore the 
magneto- Josephson effect in junctions based on spin-orbit coupled quantum wires. In contrast to the 
topological insulator case, the periodicities of the magneto- Josephson effect no longer follow from 
an exact superconductor-magnetism duality of the Hamiltonian. We employ numerical calculations 
as well as analytical arguments to identify the domain configurations that display exotic Josephson 
physics for quantum-wire junctions, and elucidate the characteristic differences with the correspond- 
ing setups for topological insulators edges. To provide guidance to experiments, we also estimate 
the magnitude of the magneto-Josephson effects in realistic parameter regimes, and compare the 
Majorana-related contribution to the coexisting 27r-periodic effects emerging from non-Majorana 
states. 



I. INTRODUCTION 

Given their exotic properties and intriguing promise 
for topological quantum information processing,'!'^ Ma- 
jorana fermions have recently received much attention in 
the condensed-matter context .I^El Promising habitats of 
Majorana fermions include the v = 5/2 fractional quan- 
tum Hall state^ as well as topological insulator edgea^El 
or semiconductor quantum wire^Silfl. proximity coupled to 
s-wave superconductors. Several experiments may have 
already provided evidence for Majorana bound states in 
semiconductor quantum wires .l^iffi^ One of the most di- 
rect but also challenging experimental confirmations of 
the existence of Majorana bound states would be based 
on the periodicity of the Josephson effect. For junctions 
of topological superconductors, the Josephson effect is 
predicted to be Att periodic in the phase difference of the 
order parame ter, in sharp contrast to the conventional 
27r periodicity^EZESI (ggg also Refs. [T9H25l for more recent 
works). 

Recently, it was noticed that a topological-insulator 
edge, proximity coupled to an s-wave superconduc- 
tor e xhibit s an exact superconductivity-magnetism 
duality.E2i2ZI xhe duality transformation maps the phase 
of the superconducting order parameter to the direction 
of the applied magnetic field in the plane perpendicular 
to the spin-orbit field. As a consequence, the duality pre- 
dicts a magneto-Josephson effect by which a rotation of 
the magnetic field across a junction induces a Josephson 
current even in the absence of a phase gradient .'^SEIl 

Explicitly, proximity-coupled topological-insulator 
edges are described by the Bogoliubov-de Gennes 
HamiltoniarPl 

"Hti = vpT^'cT^ - fir'^ + A (cos (f) t'-^ — sin (f> t^) 

-6(7^ + 5 (cos 61 cr^ - sin 9 a^) . (1) 



Here we have employed the Nambu spinor basis ^^ — 
("0^, ^^, "01, — i/'I) and introduced Pauli matrices a'^ and 
r" that act in the spin and particle-hole sectors, respec- 
tively. The edge-state velocity is given by w, p is the mo- 
mentum, and the cr^-direction represents the spin-orbit- 
coupling axis. We allow the chemical potential /i, su- 
perconducting pairing Ae*"^, longitudinal magnetic field 
strength 6, transverse magnetic field strength B, and the 
transverse-field orientation angle 9 to vary spatially. This 
Hamiltonian takes the same form upon interchanging the 
magnetic terms {6, B, 9, tr"} with the superconducting 
terms {/x. A, </), r°}. An important aspect of this dual- 
ity is that it maps the two topologically distinct phases 
of the model into each other, mapping the ' A-phase' (oc- 
curring for A^ — 6^ > maxjS^ — /i^, 0}) into the 'B-phase' 
(occurring for B^ — fi^ > max{A^ — 6^, 0}) and vice versa. 

For topological-insulator edges, the duality immedi- 
ately allows one to derive the periodicity of the magneto- 
Josephson effect from the known periodicities of the 
Majorana Josephson effect .I^ZHSSl To this end, we con- 
sider three-leg junctions^ni with the phase arrangements 
B—A—B (with a 27r-periodic Majorana Josephson effect) 
and A — _B — A (with a 47r-periodic Majorana Josephson 
effect). The duality implies that the periodicities are 
reversed for the magneto-Josephson effect, which is 47r 
periodic in the magnetic-field orientation for _B — A — i? 
junctions but 2tt periodic for a A — J5 — A setup. Strictly 
speaking, the duality also maps charge Josephson cur- 
rents into spin Josephson currents. At first sight, this 
may suggest that a change in direction of the magnetic 
field across a junction only drives a spin Josephson cur- 
rent. However, it was shown in Ref. i27ithat as a result 
of the spin-momentum locking, there is also a conven- 
tional (and experimentally more accessible) charge cur- 
rent across the junction in addition to the spin current. 

While the magneto-Josephson effect has been studied 
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FIG. 1: Upper panels: Color scale plots of the low-energy 
Bogoliubov-de Gennes spectra (only e>0)ofa_B — A — B 
junction vs. (a) the superconducting phases (j>i and <j)r and (b) 
the magnetic-field directions 9i and 9r. Lower panels: Cor- 
responding line cuts along (j)r = and Or — 0, respectively. 
Both line cuts exhibit cusps at zero energy, reflecting pro- 
tected zero-energy crossings and thus, the Josephson current 
is 47r periodic in both and 6. The parameters for the three 
segments are A;/r — 1.6, Am ~ 2, B;/,. — 2, Bm ~ 0.9, 
P'l/m/r ~ 0, and Lj = 2 (length of the junction). We also set 
9i = 9r = 7r/2, 61m = in (a) and (^i = (p,- = 7r/2, (^m = 
in (b). Note that the parameters of the three segments are 
labeled by subscripts I, m, and r. 



in some detail for topological insulator edgeSjl^StSSl much 
less is known about it for junctions based on semicon- 
ductor quantum wires. There are several reasons why 
this poses an interesting problem. Many of the ongoing 
searches for Majorana fermions are based on quantum- 
wire based structures. There are also several distinct 
differences between topological superconducting phases 
based on proximity-coupled topological insulators and 
semiconductor quantum wires. First, the kinetic energy 
of the quantum-wire Hamiltonian explicitly violates the 
duality, making the duality only of suggestive value for 
the quantum-wire situation. Second, the two topolog- 
ically distinct phases of the topological insulator effec- 
tively trade places in the quantum wire. For instance, 
a 47r-periodic Majorana Josephson effect occurs in the 
A — B ~ A arrangement in topological insulators, but in 
the B — A — B arrangement in quantum wires. 

This motivates us to explore the magneto- Josephson 
effects in semiconductor quantum wires in more detail in 
this paper. In Sec.lllj we present numerical results based 
on a recursive scattering-matrix approach and establish 
the periodicities of the magneto- Josephson effects. In 
Sec. |III[ we provide further insight into the periodici- 
ties by analytical arguments and the analysis of limiting 
cases. Finally, Sec. |IV| is concerned with numerical esti- 
mates of the magnitude of the effect and Sec. |V] collects 
our conclusions. 



II. NUMERICAL RESULTS 

We now turn to semiconductor quantum wires proxim- 
ity coupled to s-wave superconductors. The Hamiltonian 
for a clean, single-channel semiconductor quantum wire 
(QW) in the presence of a Zeeman field B, Rashba spin- 
orbit coupling M, and induced superconductivity A l^SEfll 



Hqw = ( f^ " '^ ) ^- + "^'^^^^ + ^ (^'^^+ + e"'^^-) 



A(e'"*r+-Ke-*'^T_) 



(2) 



Other than dropping the longitudinal magnetic field term 
for lack of relevance in the following, this Hamiltonian 
differs from that of the topological insulator edge in 
Eq. M by the kinetic term p^ /2ra. This term explic- 
itly breaks the duality present for the topological in- 
sulator edge and is responsible for key differences be- 
tween the topological insulator edge and the quantum 
wire. Most importantly, the phases are in some sense 
effectively reversed in the two systems. Explicitly, in 
quantum wires, the topological (or B) phase occurs for 
B^ > A^ -I- /i^, while the nontopological (or A) phase 
requires B^ < A^ + /i^. In the quantum wire model, the 
identification of topological and nontopological phases is 
unique since the A-phase is continuously connected to the 
vacuum. The corresponding identification is less defined 
for the topological insulator as the model does not con- 
nect naturally to the vacuum due to the linear spectrum. 
Indeed, the duality of the model maps the two phases 
into each other, suggesting that they are topologically 
distinct but cannot be labeled as topological and non- 
topological. However, if we take the presence or absence 
of the fractional (47r-periodic) Majorana Josephson effect 
as the defining feature of a topological superconducting 
phase, we would crudely label the A phase as topological 
and the B phase as nontopological, which just reverses 
the assignments for the quantum wire model. 

We consider quantum-wire junctions consisting of 
three segments, with phase arrangements B — A — B and 
A — _B — A. It is well established that the periodicity 
in the superconducting phase difference across the junc- 
tion is 477 in th e B — A~ B arrangement, but 2tt for the 
A — B — A setup. ^ ^'^ These periodicities are reproduced in 
our numerical calculations of the low-energy Bogoliubov- 
de Gennes spectra shown in Figs, fija) andl2ja). Here, 
we restrict ourselves to fi — for simplicity. The calcu- 
lations are based on a scattering-matrix approach which 
has been employed previously in the context of topolog- 
ical superconducting phases and Majorana fermions in 
quantum wires (see, e.g., Refs . [30l and [31]) . In short, it is 
based on concatenating small slices of quantum wire to 
obtain the scattering matrix S{e) of the entire wire. The 
spectrum can then be determined by solving the equa- 
tion det[l — S{e)] = 0. A more detailed description of the 
method can be found in Ref . [501 

Fig.[IJa) shows the low-energy spectrum oi a. B — A — B 




FIG. 2: Upper panels: Color scale plots of the low-energy 
Bogoliubov-de Gennes spectra (only e>0)ofaA — B — A 
junction as function of (a) superconducting phases and (b) 
magnetic field directions. Lower panels: Corresponding line 
cuts along 0^ = and 6,- — 0, respectively. There are no 
zero-energy crossings and the current is 27r periodic. The 
parameters for the three segments are Aj/^ = 2, Am = 1.6, 
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junction as function of the phases of the superconduct- 
ing order parameters of the outer segments. The left and 
right segments are chosen much longer than the coher- 
ence length so that the Majorana bound states at the 
outer ends do not couple to the Majoranas at the junc- 
tion and can be safely ignored. The low-energy spectrum 
shows protected zero-energy crossings, which makes the 
current An periodic as a function of 4>. The corresponding 
spectrum of a A — _B — A junction is shown in Fig. l2ja). 
In contrast to the B~A — B case, the current is always 27r 
periodic. For both types of junctions, there are two Ma- 
jorana bound states at the interfaces between the B and 
A regions. However, in the latter case the hybridization 
of the Majoranas does not generate a protected crossing 
at zero. 

Representative results for the dependence of the low- 
energy Bogoliubov-de Gennes spectra of the junctions on 
the directions of the magnetic fields are shown in Figs, 
[ijb) and 2]jb). We find that also the dependence on 9 is 
47r periodic for B — A — B junctions and 2tt periodic for 
A — i? — A junctions, with the spectra exhibiting pro- 
tected zero-energy crossings in the first case, but not in 
the second. We summarize the periodicities for the two 
types of junctions in quantum wires in Table llj Remark- 
ably, the magneto- Josephson effect has the same period- 
icity for the quantum wire and the topological-insulator 
edge. This is in stark contrast with the ordinary Joseph- 
son current which has different periodicities in the two 
models reflecting the reversed roles of topological and 
nontopological phases. 



TABLE L Periodicities of the Josephson energy as function 
of the phase difference of the superconducting gap and the 
relative magnetic-field orientation d. We list results for 
B — A — B and A — B — A junctions realized in quantum 
wires (QW) and topological insulator edges (TI). The latter 
results are taken from Ref. 1271 



B- A- B 



B-A 



periodicity for QW 47r 47r 27r 2n 

periodicity for TI edge 27r 47r 47r 2n 



III. LIMITING CASES AND ANALYTICAL 
CONSIDERATIONS 

To gain more insight into the periodicities of the 
Josephson effects summarized in Table IT] and their re- 
lations, we now combine analytical arguments and an 
analysis of limiting cases. First, we use analytical argu- 
ments to derive the Josephson periodicities for quantum 
wires which are based on the well-established result that 
the dependence on the superconducting phase is 47r pe- 
riodic for a. B — A — B junction (fractional Josephson 
effect). This complements the arguments based on the 
magnetism-superconductivity duality for the topological 
insulator edge. 

To gain a better understanding of the similarities of 
and differences between the topological insulator and 
quantum wire cases, we then study the limit of large spin- 
orbit coupling ego = mu^ for the quantum wire model 
(or equivalently large mass to), i.e. eso ^ A S> |-B — A|. 
In this limit, there are strong similiarities between the 
low-energy spectra of the topological-insulator edge and 
the quantum wire. 



A. Analytical argument 

In this section, we derive the periodicities in both 
(f> and 9 for the quantum-wire case by analytical argu- 
ments. Our arguments assume the well-established frac- 
tional Josephson effect (i.e. a 47r-periodic (j) dependence) 
for a i? — A — _B junction and reproduce the numerically 
obtained periodicities summarized in Tablell] To this end, 
it suffices to derive the number of protected zero-energy 
crossings considering particular limiting cases of the two 
types of junctions. By adiabatic continuity, these period- 
icties must then hold for junctions of the same kind with 
arbitrary parameters. 

Consider first a quantum- wire junction in the A— _B— A 
configuration. The A phase of the quantum wire is adia- 
batically connected to the vacuum by making the chem- 
ical potential large and negative. At the same time, the 
B phase is adiabatically connected to a spinless p-wave 
superconductor by taking the limit of large Zeeman field 
BW^ Consequently, a A — i? — A junction can be adia- 



batically deformed into an essentially finite segment of 
a p-wave superconducting wire with hard- wall boundary 
conditions. In this limit, the two Majorana bound states 
localized at the domain walls hybridize and split by some 
finite energy. Clearly, the two Majorana bound states 
will penetrate only very little into the A sections of the 
junctions, and consequently, they will be only weakly de- 
pendent on (j) and 9 as long as \B\, |A| <C |/^|, where /i 
is the chemical potential in the outer A segments of the 
wire. Thus, while there will be a variation of the energy 
splitting with and 6, it will be small compared to the 
magnitude of the splitting itself. Thus, there are no zero- 
energy crossings in this case and the Josephson current 
is 27r periodic both in 9 and (f>. 

We now turn to the B — A~ B junction for which the 
Majorana energy is Att periodic both in ip and 9. The An 
periodicity as a function of rep resents the well-known 
fractional Josephson effect .'^'^^'^ In the remainder of this 
section, we demonstrate that the parities of the number 
of protected zero-energy crossings of the Majorana energy 
dispersion as a function of (f) and 9 are equal for a B — 
A — B quantum wire junction. The basic observation is 
that we can again consider the limit in which the middle 
A section has a large and negative /i. In this insulating 
limit, the gap does not close when we take B and A 
equal to zero. In effect, we can thus replace the B — A — B 
junction by a i3—/—_B junction, where the middle section 
is a conventional insulator (/). 

We start by considering a, B — I interface between a 
B dominated phase with (j),9 — and a normal insula- 
tor with B, A = and /i < 0. This interface harbors 
one zero-energy Majorana bound state with wavefunc- 
tion ip. We can tune the left region to the phase 4> and 
the angle 9 by performing the unitary transformation 
U{(j>,9) = exp{i(l)Tz/2 + i9az/'2) on the Majorana wave- 
function, i.e., ip{4),9) — U{(j),9)ip. It is crucial for our 
argument that we can effect the variation of and 9 in 
the left region by a global transformation U{4>,9), which 
is possible because the rotation of B and A does not 
affect the normal insulator on the right. We note that 
{/(27r, 0) = [/(O, 27r) = —1, which guarantees that the 
Majorana wavefunction evolves to the same final state, 
when either (j) oi 9 advance by 27r. 

We now consider weak coupling of two such interfaces 
in a B — I — B junction. This coupling leads to a symmet- 
ric splitting of the two Majorana states about zero energy. 
When the coupling across the junction is sufRciently 
weak, we can obtain this subgap spectrum emerging 
from the Majorana modes localized at the junction accu- 
rately from first-order perturbation theory. Starting at 
(j) = 9 = and tuning either or to 27r the initial wave- 
function evolves to the same final state. Consequently, 
the initial and final subgap-energy spectra emerging from 
the hybridized Majorana modes will be identical for both 
processes. We know from the fractional Josephson effect 
that the positive-energy excitation at cf) — becomes neg- 
ative (and vice versa) when advances by 27r, and hence 
the associated Bogoliubov-de Gennes eigenenergy must 




FIG. 3: Bulk spectrum of the quantum-wire Hamiltonian S 
in the limit of strong spin-orbit coupling. 



cross zero energy an odd number of times in the process. 
Given that U{2t:,0) = U{0,2n) = -1, this immediately 
implies that the positive- and negative-energy excitations 
also exhibit an odd number of zero-energy crossings when 
tuning 9 from to 27r instead, which proves the An peri- 
odicity as a function of 9. 

It is worth noting that this argument fails for the TI 
edge model (fTl) , as it should according to Table IT] The 
reason is that irrespective of /i, the corresponding spec- 
trum is never gapped when setting S, A = 0. 



B. 



Strong spin-orbit coupling (eso ^ B > A) 



The arguments in the previous subsection explain the 
periodicities of the magneto- Josephson effects for semi- 
conductor quantum wires. When combined with the du- 
ality arguments for topological insulator edges, this ex- 
plains the full set of periodicities collected in Table [l] 
However, it does not fully explain the distinctly different 
periodicities of these two systems. We now address this 
issue by considering the limit of large spin-orbit energy. 

When the spin-orbit energy is much larger than the 
Zeeman energy, the bulk spectrum of Eq. ([2| depicted in 
Fig. |3]has three minima located at p = and p = ip^?, 
where pp = 2mu when /z = 0. Since pp is large in 
the limit of strong spin-orbit coupling, the subspaces at 
p — and at p ^ ±Pf effectively decouple for suffi- 
ciently smooth domain walls and the low-energy spec- 
trum can be understood as arising from a superposition 
of two subspectra.'^S Near p — 0, the Hamiltonian ([2| can 
be linearized and reduces to the Hamiltonian of the topo- 
logical insulator edge (i.e., Eq. (fTl) with fi = b = 0). Near 
p = ±pp, the Hamiltonian can be linearized, as well, and 
reduces to that of a spinless p-wave superconductor (cf. 
App.[A|). This describes a topological superconductor by 
itself. Since the topologically distinct phases are labelled 
by a Z2 index, this provides an explanation for the ef- 
fective reversal of phases between the quantum- wire and 
the topological-insulator Hamiltonian. 

The high-momentum subspace near p — ±pp has a 
gap of size A. In contrast, the low-momentum subspace 



0.08 



0.06 




0.04 



0.02 



0.5 



0.4 



0.3 



0.2 



0.1 




m=1.5 

m=5 

m^50 



0.02 




2n 



37t 



4ji 



FIG. 4: Low-energy spectrum oi a. B — A — B junction as func- 
tion of (l)r for different masses m (and hence spin-orbit energies 
eso = mv?). For large eso, the dispersion of the hybridized 
Majorana bound states becomes 2n periodic while additional 
Andreev bound states cross zero energy at cf>r — tt, Stt. The 
avoided crossings between the Andreev and Majorana bound 
state excitations vanish in the limit eso S> B. The inset shows 
that the zero-energy crossing at = tt persists for all values 
of m. Parameters: 5;/^ = Am — 2, Aj/^ — Bm = 1, u = 1, 

fJ-l/rn/r = 0, Lj = 2, and 6l/m/r = 4>l/m ~ 0. 



FIG. 5: Low-energy spectrum of a A — B — A junction as 
function of 0,. for different masses m (and hence spin-orbit 
energies eso ~ mu^). As indicated by Table hi the spectrum 
approaches 47r periodicity for eso -^ oo. For strong spin-orbit 
coupling, Andreev bound states appear. Their energy has a 
phase dependence e^ ~ Acos(0/2). Inset: Close-up near 
(f} — TT showing the avoided crossing at e = 0. Parameters: 
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near p = has a gap equal to |i? — A|, which is con- 
trolled by the competition of Zeeman and pairing ener- 
gies and which is much smaller when the system is close 
to the topological phase transition, |i3 — A| ^ A. Ze- 
roes of the gap in the low-momentum subspace trigger 
the topological phase transition and thus, the Majorana 
bound states, localized at domain walls between B- and 
A-dominated regions, predominantly reside in this sub- 
space. Consequently, in this subspace, the relevant sub- 
gap spectrum of a short junction is determined by the 
hybridization of the Majorana bound states and the pe- 
riodicities as a function of (/) and 9 are those for the 
topological-insulator edge. This seems consistent with 
Table [l] for the dependences on 6, but not for those on (p. 

To understand the full set of periodicities in Table|lj we 
thus need to also consider the high-momentum subspace 
at ±pf ■ In this subspace where the Hamiltonian reduces 
to that of a spinless p-wave superconductor, the effective 
spin-orbit field is large and hence, the magnetic field is 
only a small perturbation. The corresponding spectrum 
should thus depend only weakly on 9. At the same time, 
variations in cf) can result in a considerable Josephson 
current. In fact, as the high- momentum subspace by it- 
self constitutes a model of a topological superconductor, 
the (j) dependence of the corresponding Bogoliubov-de 
Gcnnes spectrum exhibits protected zero-energy cross- 
ings. 

We now use these insights to understand the similari- 
ties and differences of the Josephson periodicities in topo- 
logical insulators and quantum wires in more detail. First 
consider a quantum wire in a S — A — _B configuration. 



Such junctions exhibit a 47r-periodic Josephson current in 
the superconducting phase, with a protected zero-energy 
crossing of the Bogoliubov-de Gennes spectra. This con- 
trasts with the 27r periodicity for the same junction made 
of topological-insulator edges. To understand this dif- 
ference in periodicity. Fig. HI shows how the low-energy 
spectrum changes with increasing spin-orbit energy. As 
expected based on the general arguments above, the spec- 
trum develops two distinct types of subgap states as the 
spin-orbit energy increases, eso ^ B (seen most clearly 
in the traces for m = 1000 in Fig. [4]). The first type of 
state has an approximately sinusoidal dependence, an 
offset from zero energy, and 2tt periodicity. This state can 
be identified with the hybridized Majorana bound states 
in the low-momentum subspace. The second type of state 
crosses zero energy at cj) = tt, Stt with a ztEo cos{(f>/2) dis- 
persion, where Eq is of the order of A. This excitation 
corresponds to an Andreev bound state at p = ±Pf- 
As seen in Fig. |4J there is an avoided crossing between 
these states which disappears as the spin-orbit energy 
and, with it, the momentum mismatch diverge. 

This now allows one to understand the periodicities 
of Table |T] for the case oi B — A — B junctions. In the 
quantum wire, only the low-momentum subspace has an 
interesting 9 dependence. Thus, the 9 dependence re- 
mains the same between quantum wires and topological- 
insulator edges. At the same time, both subspaces con- 
tribute to the dependence on </>. Indeed, the above con- 
siderations show that the protected zero-energy crossing 
in the quantum-wire spectrum is associated with states 
which converge entirely on the high-momentum subspace 
as the spin-orbit energy increases. These states do not 



exist for the topological-insulator edge whose (/) depen- 
dence is thus 27r periodic. 

In a A — B — A junction, the change of periodicities is 
opposite. While the quantum wire is 27r periodic in 0, the 
topological-insulator edge is An periodic. The evolution 
of the low-energy spectrum for the quantum wire with 
increasing spin-orbit energy is shown in Fig. [5] The si- 
nusoidal 27r-periodic dependence of the Majorana states 
present for eso — B becomes a ±cos{4>/2) dispersion 
with avoided crossing at n and Stt for eso ^ B. In 
the limit of large spin-orbit energy, these states reside in 
the low-momentum subspace and reflect that the topo- 
logical insulator model displays a topological Josephson 
effect in this subspace. Similarly, there are also Andreev 
states in the high-momentum subspace, similar to the 
ones in B — A — B junctions with the same zLEq cos{ip/2) 
dispersion. At large but finite values of the spin-orbit 
energy, the levels in the low- and high-momentum sub- 
spaces mix, resulting in avoided crossings at = tt, Stt 
and a 27r-periodic spectrum. The avoided crossings close 
as the spin-orbit energy diverges, explaining the differ- 
ence in (j) periodicities of quantum wire and topological 
insulator. Finally, the absence of change in the 6 depen- 
dence between quantum wire and topological insulator 
has the same explanation as for B — A — B junctions. 



IV. MAGNITUDE OF THE 
MAGNETO-JOSEPHSON EFFECT 

In experiments aimed at detecting the 47r-periodic 
Josephson effect, a 27r-periodic background current origi- 
nating from the conventional Josephson effect of the con- 
tinuum states may mask the signature of the unconven- 
tional Josephson current. In the following, we provide 
quantitative estimates for the 47r- and 27r-periodic con- 
tributions to the current and show that the magneto- 
Josephson effect may be favorable over the conventional 
Josephson effect with regard to the relative magnitude of 
2TT- and 47r-periodic currents. 

In order to obtain quantitative estimates, we con- 
sider a junction with a conventional insulating bar- 
rier between two semi-infinite quantum wires in the B- 
dominated phase. In the barrier, we set _B = A = 
and /i = —Vo < 0, so that there are no unconventional 
47r-periodic Josephson currents originating from splitting 
Cooper pairs in the barrier into the two topological su- 
perconducting phases on the left and right. ^^ Thus, the 
Josephson currents in this setup are only due to the phase 
difference (f) = (pi ~ (pr or the difference 6 ^ 9i — 9r 
in magnetic-field orientations of the left and right bank. 
The total energy E (and hence the Josephson current) in- 
cludes contributions from the above-gap continuum and 
the Andreev bound states with a 27r-periodic dispersion 
(jointly refered to below as continuum contribution for 
brevity) as well as Majorana bound states whose energy 
is An periodic. 

For junctions with a low transmission probability D <C 
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FIG. 6: Maximum Majorana Josephson current Jm for an 
SIS junction normalized by the normal-state transmission co- 
efficient -\/l) of the junction as a function B and A. We set 
fi = in the superconductor, which thus supports a topo- 
logical phase for B > A. Inset: cuts along the dashed 
and solid lines in the color scale plot. The parameters are 
eso = O.OSmeV, m = 0.015me, Vo = 250meV, Lj = 3.2nm. 



1, we find the energy of the Majorana states to be 

^Majorana (0, &) = Em COs(0/2) COs(0/2 + ^o (</-))■ (3) 

The 9 dependence of the energy involves a phase shift, 
whereas the dependence is always symmetric with re- 
spect to (j) = (cf. the lower panels of Figs. [1] and [2]). The 
largest energy splitting is given by Em ^ vDEgi^p, where 
-Egap denotes the magnitude of the gap in the two banks 
of the junction. The size of the splitting is determined 
by the single-electron tunneling amplitude oc ^/D. 

The critical current of the junction depends on 9 with 
a maximum critical current of Jm = {e/fi)EM- In Fig. pi 
we show numerical results for Jm, normalized by vD 
to make the results insensitive to detailed properties of 
the tunnel junction, as a color scale plot. This normal- 
ized Majorana current roughly corresponds to (e/fi)£'gap. 
Thus for a fixed A, there is an optimal value of B for 
which the ratio Jm /vD is maximized. 

In Fig. [6J we compare the amplitudes of the 4tt- and 
27r-periodic currents originating from the Majorana and 
continuum states, respectively. We distinguish between 
the modulation of the continuum current with (p and 9. 
According to our numerical results, the phase dependence 
of the continuum energy can be well described by 

Econt{(f>, 0) ^Eo cos((/.) {a -f (1 - a) cos [9 + 9i{(j))]} 

-I- const., (4) 

where a G [0, 1] is a parameter-dependent constant. 
Hence for a > 0, the amplitude of the oscillations as 
a function of 9 is smaller than the amplitude of the (p- 
oscillations. We numerically calculate the largest ampli- 
tudes of the Josephson current oscillations as a function 
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FIG. 7: Numerical results for the Majorana {Jm) and con- 
tinuum (J0, Jg) current oscillation amplitudes as a function 
of B. For presentation the continuum current has been multi- 
plied by 30. The calculations have been performed with the 
same parameters as in Fig. [6] with A = 5fieY. Transmission 
probabilities range from 0.1% to 0.4%. The vertical dashed 
line denotes the critical magnetic field. 



of d and 



J0 = eEo/h 
Jg^{eEo/h){l-a) 



(5) 



and plot them in Fig. [7] as a function of B along with 
the amplitude of the 47r-periodic Majorana current Jm- 
The latter is much larger than the continuum contri- 
bution for a large range of parameters. For a p-wave- 
superconductor junction the 27r-periodic current involves 
tunneling of Cooper pairs with amplitude ex Dj^in con- 
trast to single-electron tunneling ex vl? responsible for 
the 47r-periodic current. Hence in the large B limit, we 
expect Jm to exceed the 27r-periodic current by a factor 
of 1/vD, which is ~ 20 for the parameters used in Fig.M 

Only very close to the phase transition, when \B — 
A| ^ A, the continuum current can exceed the Majo- 
rana contribution. This is consistent with numerical es- 
timates for continuum and Majorana Josephson currents 
for the topological insulator edge (IT]) in Ref. [201 which 
corresponds to the limit |B — A| <C A for the quantum 
wire model pi. 

Comparing the Josephson and magneto- Josephson ef- 
fects, we find in accordance with Eq. (|5|) that J^ is larger 
than Jg, in particular in the regime of small B. On the 
other hand, the amplitude Jm of the Majorana current 
oscillation is the same for (f) and 6. Thus, near the crit- 
ical magnetic field, the 47r-periodic magneto- Josephson 
current appears on top of a constant current background 
with a small 27r-periodic modulation from the continuum 
states (see Fig. Isk) . This is favorable in experiments to 
discriminate the 47r-periodic Majorana current from the 
conventional Josephson current of the continuum, e.g., in 
the Shapiro-step-like pattern due to the interference of a 
rotating magnetic field and an ac voltage as described in 
Ref. [571 
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FIG. 8: Numerical results for the continuum contribution to 
the Josephson current as a function of </> (red crosses) and 
Q (green dots). The solid lines are fits to the (/i-derivative 
of Eq. (4), d^Econt{4',S)- All parameters are the same as in 
Fig. [T] The fixed angle is chosen such that the current is 
maximized. In this way, the oscillation amplitudes can be 
used to obtain J^ and Je plotted in Fig. [7] The dependence 
is simply 7(0) — J4, sm{4>). In contrast j(d) has a large offset 
and only a weak 9 dependence for B ~ A [panel (a)]. This 
corresponds to a ~ 1. For B 3> A [panel (b)], the amplitude 
of Jg approaches that of J^ and a decreases. 



V. CONCLUSIONS 

The dedicated experiment^iil^i^ to detect Majorana 
bound states in spin-orbit coupled quantum wires raise 
the question how the exotic signatures of Majorana zero 
modes manifest themselves in such wires. In this work, 
we explored the magneto- Josephson effect which comple- 
ments the remarkable Josephson physics that Majorana 
bound states entail. 

Our principal goal was to determine the periodici- 
ties of the magneto- Josephson effect for the various do- 
main configuration of the Josephson junction. For junc- 
tions made of topological insulator edge states, the 47r- 
periodic effects emerge in mutually exclusive configu- 
rations: the conventional Josephson effect (involving a 
phase difference of the superconducting order parameter 
across the junction) appears in the A — B — A domain se- 
quence, while the magneto- Josephson effect requires the 
complementary structure, B — A — B. This indicates 
that in a sense, both domain types are topological for 
a topological-insulator edge. In contrast, for spin-orbit 
coupled wires, we found for both types of Josephson ef- 
fects that a 47r periodicity requires the B — A — B con- 
figuration. 

While we invoked both analytical arguments and nu- 
merical analysis to establish this result, it has a natural 
interpretation. We expect that 47r periodicity with a pa- 
rameter which is normally defined between and 27r (up 
to trivial shifts) can only emerge if the parameter per- 
tains to a topological phase. In a spin-orbit coupled wire, 
there is no ambiguity as to which phase is topological. A 
A-dominated phase is continuously connected to the vac- 
uum by taking the limit of a large and negative chemical 
potential. The i?-dominated phase, on the other hand, is 
a topological phase continuously connected to a spinless 
p-wave superconductor. In Sec. [Ill A] this argumentation 
is made explicit using analytical arguments. 

Despite these characteristic differences of the magneto- 



Josephson (as well as Josephson) periodicities between 
topological insulator edge and semiconductor quantum 
wire, both models can be connected explicitly in the limit 
of large spin-orbit coupling. We exploited this connec- 
tion in Sec. lIIIBJ to understand the relation between the 
Josephson periodicities of the two models explicitly. 

With a view towards experiments on Majorana Joseph- 
son phenomena, we also computed the 47r-periodic 
magneto- Josephson effect for typical parameters, and 
compared its magnitude to that of the more conven- 
tional 27r-periodic background. The above-gap contin- 
uum of states in the wire contributes to both the phase- 
controlled and the magneto- Josephson effect. In the low- 
transmission regime {D <^ 1), we found that both Att- 
periodic Josephson effects yield currents of the order of 
j^\/DEgap, with Egap the gap in the two banks of the 
Josephson junction. In contrast, the conventional effects 
are suppressed by an additional factor of vZ). Experi- 
mentally, these effects may be easier to observe for larger 
transmission probabilities D. In this regime, however, 
the exotic and conventional current contributions are of 
the same magnitude, although the 47r periodicity is rel- 
atively more pronounced for the magneto- Josephson ef- 
fect. 

Josephson-related phenomena in spin-orbit coupled 
wires can become more complex when considering, e.g., 
ac modulations and Shapiro steps. These may require 
one to take accurate account of the complicated spectrum 
of the quantum wires. Moreover, when the energy gap in 
the middle domain is not too large, additional Andreev 
bound states could be present, which contribute to the 
Josephson effect. In this work we refrain from discussing 
these topics as well as more complicated setups such as 
three-leg Josephson effects to keep the presentation con- 
cise. Nonetheless, these aspects may prove important 
(and maybe even beneficial) in experiments, and present 
interesting avenues for future research. 
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Appendix A: Quantum wire in the limit of large 
spin-orbit coupling 



In Sec. |IIIB[ we relied strongly on the statement that 
in the limit of strong spin-orbit coupling, the quantum- 
wire model (pi) reduces at low energies to a combi- 
nation of the topological-insulator-like low-momentum 
subspace and a spinless-p- wave-superconductor-like high- 
momentum subspace. In this appendix, we provide an 
explicit justification for this statement. 

The statement is evident for the low-momentum sub- 
space, so we will not consider it further. For /i = 0, the 
Fermi points are located at p — (low momentum) as 
well as p — ±pf with pp — 2mu (high momentum) . As 
the spin-orbit coupling (or equivalently, m) increases, pp 
becomes large and so does the effective spin-orbit field in 
the high-momentum subspace. Thus, in this limit, we can 
first diagonalize the Hamiltonian in the absence of the in- 
duced pairing A and then treat the latter perturbatively. 
To do so, we perform the unitary transformation 



U = exp{iaayT2/2} exp{idaz/2} 



(Al) 



on the wire Hamiltonian (pi) with A = 0. If we choose a 
such that tana = B/up, the rotated Hamiltonian takes 
the form 



Ho 



2m 



+ y/{wpfnP(T, 



(A2) 



The low-energy subspace at p — ±pF is formed by the 
bands for which a^ takes the value —1. We now reintro- 
duce the pairing term At^ and apply the transformation 
U to it. The projection of "Ho onto the lower bands yields 



^cff — 



P 
2m 



- ^/{upY + bAtz + 



upA 



B2 



(A3) 



The condition £50 ^ A guarantees that we can neglect 
the coupling to high-energy degrees of freedom near ±pf- 
Linearizing around the Fermi momenta and using \up\ ^ 
^so ^ B the effective Hamiltonian takes the form 

"Hcff = u {\p\ - pf) t^ -I- sign(p)ATa;. (A4) 

This describes a spinless p-wave superconductor. 
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